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1.0 INTRODUCTION 

This is the final engineering report on Task D of NASA Contract No. 
NAS6-2307. This task comprised a study of experiment requirements, of 
technical characteristics, and of the GEOS-C radar altimeter related analyses. 
In addition to the work reported herein, a study of engineering test data re- 
quirements was also conducted; results of this . activity were documented in a 
report distributed in December, 1973. 

Chapter II of this report contains statistical analyses related to deter- 
mination of wave height resolution achievable as a function of system character- 
istics and averaging period. An equally important topic of this chapter is the 
desirability of using computer procedures to compensate for altitude tracker 
time-jitter. 

Chapter III examines data processing considerations for the GEOS-C 
system. An extensive analysis of the spatial filter effect is given and 
results of a computation of geoidal power spectral density, based on Skylab 
altimeter data, is displayed and interpreted in terms of projected GEOS-C 
random errors. This information is then used in deriving minimum-mean- square 
filter procedures for both geoid undulation and slope data. 

Chapter IV examines the characteristics of mean received waveforms as a 
function of off-nadir angle. This information is then used to obtain tracker 
bias as a function of sea state and pointing angle. The angle estimation 
process proposed by the GEOS-C hardware contractor (General Electric) is 
also investigated from a standpoint of achievable angular resolution. 
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2.0 ANALYSIS OF GEOS-C WAVEHEIGHT RESOLUTION 

2 . 1 Summary of Results and Recommendations 

The main body of this section, is devoted to an analysis of waveheight 
resolution, in significant waveheight units, achievable with the GEOS-C 
altimeter, and the effect of uncorrected altitude tracker time-jitter on 
resolution. The final results of this analysis are shown in graphic form in 
Figures 2.1 and 2.2. Results for pulse lengths of 10 and 12.5 ns are given 
because currently available GEOS-C test data shows the pre-detection pulse 
width to be ~ 13 ns at the 6 dB points. Since a square-law detector is used, 
this implies a post— detection pulse width of 13 ns at the 12 dB points and 
we estimate that the video filter would increase this pulse width to ~ 10 ns 
at the 3 dB level. Figure 2.1 shows that for the 10 ns pulse length with 
expected t im e jitter (*“3 ns) , 80 percent of the observations will be within 
±25 percent of the correct value for significant wave heights of ^ 3 

meters. Below this H-^ value the resolution degrades rapidly. For a time- 
jitter of 6 ns, the corresponding H-^y^ value is - 3.5 meters. To allow a 
comparison with other results, this figure also shows a "one- sigma" confi- 
dence level result; here the 25% resolution value intersects the curve at 
lij. a 2.4 meters. We do not recommend use of this bound. 

Figure 2.2 shows similar results for a 12.5 ns pulse width. Here the 6 ns 
jitter case ±25% resolution point is seen to occur at ~ 4.25 meters. For 
an rms tracking jitter of 12 ns (which is considered to be in the category of 
substandard data) the ±25% resolution point is “'5.15 meters. Note that if the 
12 ns rms jitter is assumed to be correctable (at H^y^ - 3.6 meters) tracking 
jitter correction is equivalent to an extension of the averaging period by a 
factor of 2.5 (waveheight resolution is proportional to the square root of 
the averaging period) . 

Both Figures 2.1 and 2.2 are based on an averaging period of 20 seconds. 
This value was used since we feel it represents an approximate upper bound 
on averaging period, unless a. priori information is available regarding ocean 
surface homogeneity. .Since high seas are of quite low probability on a global 
scale, it would be highly desirable to be able to collect experimental wave- 
height data as it exists, in the planned GEOS-C calibration area. For an 
experiment that requires high seas (e.g., the North Atlantic in the winter 
months), the acquisition of aircraft sea-truth data is a formidable task. 






Waveheight resolution (percentage of H 
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To provide some insight into the problem of extending the data averaging 
period, two randomly selected sea-truth charts are shown in Eigures 2.3 and 
2.4. These were supplied by NOAA during the Skylab program’s SL-2 and SL-3 
missions. Figure 2.3 shows a 3 meter isocontour crossing for pass 9 

and 2 meter isocontour crossings for pass 8 - for the wind-driven (non-swell) 
data shown. The scaled time periods for the ground track to traverse these 
areas are 10 and 18.6 seconds. In Figure 2.4 the subsatellite path crosses 
two 1-meter contours in approximately 20 seconds. 

Based on the above discussed results, the ~3 ns rms tracking jitter 
level anticipated for the GEOS-C altimeter is seen to produce a relatively 
small degradation in wave height resolution. Therefore, under assumed nominal 
altimeter performance we do not recommend corrections in the waveform averaging 
process to compensate for the altitude tracker time-history. If, after launch, 
for some reason the altitude data quality should seriously degrade (>12 ns 
rms noise) the question of time-jitter correction should be re-examined. The 
need for corrections in the waveform averaging process for deterministic 
orbital effects is an unexamined question. For an assumed altitude rate of 
change of 50 meters /second, due to orbital eccentricity, and a tracking loop 
correlation interval of ~50 milliseconds, the altitude change during the 
correlation period is 16 ns. Hardware test data should be available in the 
near future to permit an assessment of the effect of such altitude rates on 
tracker variance. 

> 

2.2 Analysis of Tracking Time- Jitter Effects and Waveheight Resolution* 

Our main purpose in this section is to anlayze the effect of altitude 
tracker— induced time-jitter in the sample— and— hold (S & H) data on waveform 
averaging and waveheight resolution. Because of the small tracker (and thus 

i — t* * 

S (t H gate position) standard deviation (~ 3 ns) relative to the "flat sea" 
rise time, we desire to examine the effect of neglecting tracker time-jitter 
in the reconstruction process. Figure 2.5 shows presently planned waveform 
data processing procedures and the overriding importance of this effect. 

In the following analysis we examine the uncertainty in the estimate of 
■the received waveform as a function of averaging period. Since a single 
received waveshape is an ensemble member of a random process, the greater the 
number of waveforms averaged the lower the uncertainty in mean value and thus 


*A major part of this section was contributed by A. C. Nelson. 
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32/FRAME ALTITUDE AVERAGES 
320/FRAME TRACKER SERVO ERROR 
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SMOOTHED ALTITUDES, 4/FRAME 
TRACKER VARIANCE, 4/ FRAME 
16 S&H FRAME AVERAGES 
16 S&H VARIANCE ESTIMATES/FRAME 
NOISE, ATTITUDE/SPECULAR GATES 
GMT TIME TAGGING 
SELECTED HOUSEKEEPING DATA 
ETC. 


TAPE 3 CONTAINS: 

SMOOTHED ALTITUDES, 4/FRAME 
TRACKER VARIANCE, 4/FRAME * 
TIME-REALIGNED S&H FRAME 
AVERAGES AND VARIANCES 
NOISE .ATTITUDE/SPECULAR GATES 
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SELECTED HOUSEKEEPING DATA 
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03 


Figure 2.5 


Geos-C Waveform Processing for Mode 2 (first of the two high-data-rate 
telemetry modes). 1 Frame = 3.277 seconds. 
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the greater the waveheight resolution. The upper limit on averaging period 
is determined by ocean surface homogeneity; typical data relating to waveheight 
spatial variability was given in Figures 2.3 and 2.4. Receiver noise can be 
shown to produce little effect on waveheight resolution for pre-detection 
signal— to-noise ratios ^ 10 dB, and this effect will not be considered. 

In the following we will assume that the point-target pulse shape and 
the waveheight probability distribution are both Gaussian functions. This 
will give rise to an average waveform which is a cumulative Gaussian, as shown 
in Figure 2.6, and will permit closed form analysis of several important 
parameters. Since the received waveform corresponds to the integral of a 
Gaussian which in turn is the convolution of pulse shape and waveheight Gaussian 
functions, the composite variance is: 


+ a 


2 

t 


or from Figure 2.6 the composite standard deviation in nanoseconds is 

T = f (.508 H„ ,.) 2 + (.602 T ) 2 
s 7 x/3 p 

where H- is a significant waveheight in feet and T is the 3 dB video pulse 
1/ 3 F 

width in nanoseconds. Therefore, a ramp period, T r (the dotted curve in 
Figure 2.6) may be defined as: 

T = 2.58 V.258 H n , 2 + .362 T 2 
r T 1/3 P 


= 1.3iVh 1/3 2 4- 1.4 T 2 . 

This represents the time expanse of. -.the ramp period as defined by the projection 
of the mid-point slope as shown in Figure 2.6. On this basis the ramp period 
is equivalent to 1.55 T p for the case of = 0. Other rise time definitions 

could have been used. For example, rise time of pulse-like waveforms is 
sometimes taken as the 7 and 93 percent amplitude points. The definition 
used here is convenient because of its relationship to the mid-point slope. 

First we consider the statistics of tracking jitter as these are affected 
by quantization and correlation properties of the altimeter system. In the 
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Radar Pulse Shape Wave Helght Dlstrlbutlon 



2 m 

°t = 4 (H 1/3) = * 508 H 1/3 

(H 1/3 In ft.) 



for T = 3db video pulse width in ns 


Figure 2.6 Waveform Relationships 
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Skylab altimeter, altitude tracker standard deviation and quantizer step- 

size were comparable in magnitude; therefore, on the average, the output 
signal constituted roughly a 4 level discrete random process. With the 
much smaller step-size used in the GEOS-C system, the range tracker signal 
more nearly resembles a continuous- time process. 

Statistics of the mean-wavef orms derive both from the chi-squared dis- 
tribution of the basic received waveforms and from the superimposed time 
jitter of the tracking process. Since 100 statistically independent waveforms 
are received per second and the tracker jitter decorrelates in about 10 wave- 
forms, the central question is the time-wise behavior of the combined process. 

Starting with the normalized cumulative Gaussian to describe the mean 


waveform, 


<t-y)/T s 

g(t) = /e* p {- ^ <*x. 


the probability that the jitter is kx (where T = 5/32 ns is the GEOS-C 
quantization value) is given by the discrete Gaussian density function. 


= k - ±1, ±2, ±3- 


2a 


*>0 


k ■ °> 


where y is the mean on the time scale corresponding to the normalized mean 
voltage equal to 0.5, and a is the standard deviation of the jitter process. 
We first consider the idealized case in which the tracker signal only 


executes step changes every 10 pulses and then examine the validity of this 
simple model. Assume that N fc = 10 Independent samples of the tracker signal 
and N rf = 100 independent waveform samples are available per second; then 
it is desired to relate the variance of the average of N waveforms to the 
parameters T g , and ( assumed to be 3 ns), in order to provide some 

knowledge of sea state' resolution. The time varying received signal has a 
<X 2 (chi-square) distribution with 2 degrees of freedom, with mean and standard 


deviation both equal to g(t + kx) where kx is the tracking jitter. With no 
jitter the average of N = 100 independent or uncorrelated signals would have 
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2 

a variance equal to g (t)/100. With tracking jitter the average of 100 such 

waveforms would have a mean and variance as determined below. 

The following derivation is limited to the determination of- the mean 

2 

and variance because the average of fifteen X variables can be assumed to 

2 

have approximately a Gaussian distribution. Typically, the X distribution 

with 30 degrees of freedom is approximated by the Gaussian distribution. 

2 

This would occur on averaging 15 X variables each with two (2) degrees of 

2 

freedom because of the additive property of the X distribution. The 
determination of the mean is not affected by the correlation between suc- 
cessive values of the tracking jitter and the results given below are general 
in this respect. On the other hand, the correlation of the tracking jitter 
does affect the variance but this is minor when the sample size is large. 

The effect is small because: (1) jitter variance is small compared to the 

composite variance, and (2) relative length of the averaging period is 
large, (say, greater than 1 second) as compared to the less than 0.1 second 
period over which the correlation is positive. Utilizing these remarks, the 
following derivations are performed. 


2.2.1 Computation of the Mean With Jitter Present 

Let the average of r = 10 waveforms (for 0.1 sec.) be denoted by v(t) 
and for N waveforms (or N/100 seconds) by v(t). The jitter value is fixed 
for r waveforms and then shifts to another value independently selected from 
the statistical distribution. This is an approximation to the actual waveform 
process with autocorrelation function given by Table 3.1, of [1] which shows 
that on the average, the actual waveforms shift and decorrelate in~0.1 seconds 
The average of N waveforms will be given by: 


where 


= 2 m^v^/N, n = 2 


N/r, 



2 v.(t+kT)/rn, 

j=l j 


E{^} = E n E- (v | n) 


= E {2 rn, g(t+kt)/N} 
. n k. 
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= Z rnp k g(t+kx)/N = Z P k 8 k > g k = g(t+kx). 


If the jitter variable is denoted -by and the composite variable 
(waveform) by X^» the sum (for the case t — p which occurs at the midpoint 
of the ramp) 


Z p k g(t+kx) 

can be approximated very closely for small T (jitter step in ns), by the 
integral: 


/ 


Prob . { X 1 


is in the 
interval (k,k+dk) 


} • Prob . 


{ X 2 is less than } dk 
or equal to k 


and this is equivalent to: 

p{x 1 > x 2 ) 

which is equal to 0.5 for t - p because is N(0,<7.) and X 2 is N(0,T g ). 

If t B y + T say, it can s imi larly be shown that the sum Zp k g(t+kl) is 
approximately equal to 

Prob.. {X 2 <_ X x 4- T g > 

where X.. and X„ are distributed as above and thus X~-X, -T is N(-T , { a?+T 2 ), 
1 * 22 ^ _L & 

As an example, for T = 8.8, O, = 3, T +a. - (9.3) 

s j b J 


Prob. 



2 

s 


0.9462 


" 


0.8279 


Thus the mean is biased slightly downward at t = P+T ; without jitter it is 

u 

0.8413. Figure 2.7 shows the mean w^yeform distortion due to uncorrected 
tracking jitter.' 

Note that this result can be generalized to any multiple m of T g to 

yield 

£ p k g(t+kr) 

i ^ 

* Prob. (X„ < X.+mT } 
f Z — is 


This notation denotes X^ is normally distributed with mean 0 and standard 

deviation a.. 

J 



GEOS INTENSIVE MODE 


a. =' 3 ns 
J 



Jitter Free 
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given by 


where 


mT 


Prob. I u < 




2^2 
+T } 


■ l 


(0.508 H 1/3 ) 2+ (-602 T ? ) 2 


0 . = standard deviation of tracking jitter (in ns) 

3 

**1/3 ” wave ^ e:i -8ht (i n feet) 

2.2.2 Computation of the Variance With Jitter Present 

The variance of v(t) is computed by use of the formula, 

a 2 {v) = E a-{v|n} + a^E-{v[n}. 
n v 1 H v 

See reference [2] for these general formulas for the unconditional 
means and variances in terms of the conditional means and variances. The 
subscript variable indicates the one with respect to which the moment is 
being taken conditional on the remaining variables being fixed. 

E B 4{— 5^1 ■ E a 7 2 r2 V 2(t+kt)/m J 


N 

1 r 2 
m Ep A 


2 r =i i 2 S rn k v k: 

Vv {v l n} = a n E v N . 


= a. 


2 2 xti^gCt+kx) 


n 


N 


} 


.-if, 

„ 2 L 


l 8^ np k (1 “ p ^^ " 2 E s v g 0 np v p el 




k<2, ~ ~ ~ 1 

Substituting the results in equation (2) yields 

0 2 {v> = ^ H'P k 8^ 2 4 p k (1 ' p k )_ n E Wk S r 
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Using the following relationship 


Sp k S k + 2 k<£ P k P £ 8 k 8 &* 


The above becomes 


2 r=n (r-KL) „ 2 1 ^ _ s 2 

O {v} = Ip k g k - - (Sp k g k ) • 


If n * 10, r = 10, N = 100, 

2 r =, _ 11 „ 2 1 /v .2 
0 ^ 100 ~ 10 ^ P k®k * 

5 

The calculations become tedious for T - ns » as there are at least 
2 x 6 x 32/5 = 87 intervals for inclusion of 95% of the jitter values. An 
approximation was run by using t — 40/32 ns and the results extrapolated 
to t = 5/32 ns for the case in which 

T = 8.8 ns 
s 

V 3ttS 

t = u and t = ]i + T • 

In this example, at t = p . 


£ p fc g£ = 0.2643, £ p k g k = 0.5 
E{v} = 0.5. 


Therefore with jitter 

C 2 {v} *•*■‘(•064)^. 

_ 2 

Without jitter, assuming independence, c 2 {v) = (.05) . 

At t - V + T s , £ p k g^ « 0.6927, £ p^ - 0.828. 

With jitter E{v} = 0.828 

a 2 {v} = (.087) 2 . 

Without jitter, a {v} = (.*084) . 

These results are shown in Figure 2.8. 
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2,2.3 Comments 


Some comments ox remarks boncerning these results are in order. First 

suppose that the results are compared to that of averaging 100 independent 

waveforms with jitter present (no correlation). For O — 3 ns, T s =8.8 ns, 

t=p + T , and for an approximation to be subsequently discussed, 
s 


2 


0 


(1 waveform with jitter) = (corrected mean) 


+ 



_1 

Tie 


(.836)“ 


2 2 
0 (average of 100 waveforms) = (.0836) 


which may be compared with the .087 results of the previous, more exact, 
analysis. Also note that the approximation which assumes that 0 is fully 
correlated over 10 pulses 


2 

0 


.average of ^ 
''lOO waveforms'' 


(corrected mean) 

N * 
prf 


2 

0 . 

NT'-* 
T s 


_l_ 

Tie 


(* 90) 2 


£where N ^ = 100, N^ = 10J 

is seen to result in the same order of approximation. As expected, variance 

of the true process is midway -between two limiting cases of tracker signal 

2 2 

correlation. The term 0./T TTe in the above equations represents a trans- 

3 s 2 2 

formation, by Taylor series expansion, of nanoseconds to volts . Note that 

2 — 

the formula for 0 {v} is directly proportional to 1/n, i.e., the variance 
decreases as n, or the standard deviation as 1/ /H. Also for r = 1, the 
result reduces to 


a 2 {5} = f S P,ft - ^ PA) 


This result is applicable when the autocorrelation function is zero. 

The variance of v increases slowly as r increases for N fixed (=100 say) 
as tabulated below for N. = rn = 100, t = p + T . 


r 

n 

0 Z {9} 

1 

100 

.0070 

2 

50 

.0071 

4 

25 

.0072 

10 

10 

.0074 
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This shows the results to be largely insensitive to the detailed averaging 
properties of the tracking loop. 

2.2.4 Sea State Resolution 

Suppose that only significant waveheight is subject to change, then 

T « V 0. 258 H?,- + .362 T* 

s’ 1 / j V 

= ^0. 258 H^ 3 + 56.63 

for T = 12.5 ns. Assume that a. is given or has been measured to be 3 ns. 
p 1 

The following figure illustrates the effect of SWH on the total variance 
of the process. 



U li+T sl V+T s2 


The difference -of two values (for different values) is 


2 2 

TT - T 
s2 si 

' 0,258 


1/3,2 “ 1/3, 1 


2 2 

‘The difference T _ - T is estimated by the difference in the estimated 
s2 si 

values obtained from the data or by squaring the values read from the cumula- 
tive distribution curves at 0.8413 and taking their difference. The values 
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A As 

from the curve are denoted by an ^ T g ^, res P ec tively. This difference 
is dependent on the values through the following approximation: 



0.258 


(H l/3 f 2 


H l/3,1> . 



“ .0158 (H 


2 

1/3,2 


H l/3,1 } 


* 


and for 1 = 6 ft, 2 = 1, 3, 9, and 18 ft, the difference T g2 - T gl 

ranges from -*0.56 to 4.61 ns. For H^y^ i “ 9 ft and H ^y 3 2 = 1, 3, 6, and 

18 feet, T „ - T , varies from -1.29 to 3.91 ns. (See the following 
s2 si 

tabulation) . 


H l/3,1 


6 ft. 


(Approx. ) 



H l/3,1 


9 ft. 


(Approx. ) 



H l/3,2 = 1 
3 
9 

18 


-0.56 

-0.43 
0.72 
4. 61 


H l/3,2 = 1 
3 
6 

18 


-1.29 


-1.16 

-0.72 

3.91 


Now consider the precision of the estimated variances (standard 
deviations) as read from the averaged waveform. The variance of the 
ordinate v(t) was obtained earlier as a function of n (no. of 0.1 second 
intervals), p fc , and For t = y + T g , = 9 ft, the estimated variance 

was given as 0.0076 volts^* 

' 2 2 
This variance in volts can be transformed to nanoseconds by using 

f *■ ' 

the Taylor series expansion of the Gaussian $(t ; T ). (The inverse of 

os 

this was used earlier to convert jitter variance to a transformed jitter 
2 

variance in volts ). Given t = + T and 

s 


v T s ) + (t-t o ) $ <t o ; T s ) + 

0 {v} * a{t}. 

T I'frreT 


« • » 
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Thus 

a{t} = T /2ire a{v} 

S . 

o 

and using a {v} = 0.0076, the variance of a 1 second average, the following 
values are obtained for several averaging periods. 


Table 1 


a{t}, /2a{t} vs. Averaging Time 

(Seconds) 

at t = y + T , v(t) = 0 
s 

Averaging Time 


/2a(t} 

(Seconds) 

c(t} 

1 

3.17 

4.48 

10 

1.01 

1.43 

20 

0.72 

1.02 

40 

0.51 

0.72 

80 

0.36 

0.51 

120 

0.29 

0.41 

Assuming that the difference in two 

times has 

a standard deviation of 


(since the times are assumed to be equal under the hypothesis of no difference) 
then the above results need to be multiplied by /2 = 1.414, see last column 
of Table 1. Hence, for average waveform data collected over one second, the 
standard deviation of the observed difference is estimated. to be 4.48 ns. 

In comparing data for = 9 ft with that for = 1 the 


expected or average difference of T £ ” T g ^ - 1.29 ns. 


The probability that 
the hypothesis T - T will be rejected is given by the power of the test 

Si. 0 4. 

(see [3], p. 229). 

, 1.29 


1.02 


= 1.265, Power = .24 for a= .05 


i.e., for a level of significance of a = 0.05 (a 5% risk of stating that 
the H 1/v differ when in fact they do not differ), the probability that the 
H^ 3tg will be indicated as different is approximately 0.24. If X = 2.8, 
the power is 0.80, i*e. 

T s2 ” T sl = 2,8 (1,02) = 2 * 86, 


a value exceeded only by comparing ~ 9 ft vs. H^y^ “ 18 feet. 
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Similar results may be obtained using a chi-squared formulation. For 

an averaging period of 10 seconds, first assume that 10 waveforms are 

available j averaged over 1 second. From the above tabulation the standard 

deviation to be associated with a one-second average waveform is 4.48 ns. 

In order to obtain uncertainty values for an 80% confidence bound, chi- 
2 

squared (X ) values for 10 and 90% levels are used, i.e., 


2 

a MAX 



2 

2 = no 

° MIN 2 

X 90% 


9(4. 48) 2 = 9(4. 48) 2 

4.17 14.68 


a — 6.58 
MAX 


a MIN “ 3 ’ 51, 


where n = number of degrees of freedom. Using the previous equation 

t s 2 ” T sl = * 0158 ^1/3,2 ” **1/3,1^ 


with the above range of a values; as an example if one va l ue is taken 

as 8 ft, we find the other to be 16.1 ft. As another example for 1 = 

3 ft, we find 2 = 14 ft ' This means that the probability that a wave- 

height observation will lie within these values is 80 percent. 

The preceeding analysis has focused on the general statistical aspects 
of the mean waveform, variance, and waveheight estimation areas. This work 
will next be used as a framework for deriving an approximate, closed form 
solution to the waveheight resolution question. The initial problem formula- 
tion assumed Gaussian functional forms for both the system waveshape and the 
waveheight probability distribution function based on justifications given 
in [4] and [5], and because this leads to analytical tractability . Without 
the latter assumption, as must be the case with initial radar altimeter 
waveform studies, a considerably more complex deconvolution approach is 
•needed to extract waveheight information. That is, for the Gaussian 
assumption only a differentiation is required to recover the combined system 
waveshape and waveheight distribution. Except for the increase in process 
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variance which may result from differentiation, we know that the ±T g point 
on the mean waveform will directly map on the differentiated waveforms. 

For this reason the following waveheight analysis considers only the average 
received waveforms. For future reference, note that variances will propa- 
gate through differentiation (for a continuous- time analogy) as follows: 
for the signal power spectrum, which at the output of the square-law 

detector is the convolution of the IF spectrum, S(o)), the differentiated 

2 

process variance, is 

CO 

2 f 2 

a - / u) S(u)) da), 

o 

For the noise spectrum (this analysis is assumed free of' receiver noise) . 
an identical form applies. Therefore, the differentiation operation - will 
alter the signal- to-noise relationship only ta the extent that these spectra 
differ (c.f., [1] for a computation of S(u))). For essentially band limited 
processes, this factor is rather small. 

To proceed, the previous analysis showed that the waveform standard 
deviation in volts could be transformed into time uncertainty tf(t), as 

T /27re _ 0.36T 

a(t) ^ — , a(v) - 

/t /t 


t in seconds. 

Also as previously discussed, the approximate total time uncertainty 

is due to the uncorrelated variance of T along with a. ; the waveform 

s J 

statistics and the tracking jitter random sources. Denoting the combined 
one-sigma variance as cr(t), 


a 2 (t) 


(0.36 T )' 


lOt* 


Using the previously defined value for T , 


T. = .602 /712 Hr,„ 
s 1/3 


+ T 


2 

P’ 



24 


and using a differential of T g to convert the time uncertainty a(t) into 
an uncertainty gives 


AT - 
s 


• 428 H l/3 AH l/3 




712 H l/3 + T P 


2 2 

Equating (AT) with a (t) and substituting for T gives a significant wave- 

height uncertainty a of 

H l/3 


6 


256 (. 712 H^ /3 + T p ) 


2 2 2 2 2 
v + ,546a7 (. 712 H- + T ) 
3 1/3 __Jg_ 


H. 


1/3 


* 1 / 3 ^ 


This gives a first order solution to the one-sigma uncertainty in wave- 
height resolution. In accordance with previous work [6], we adopt an 80 
percent confidence bound (for a Gaussian process this bound is approximately 
1.3a). Therefore 


resolution ^ - 
(80% confidence bound) 


i 


.433 (. 712 H^ /3 + T^) 2 


H 


1/3 


+ .923a? (.712 H?. + T^) 

3 1/3 p 

/t 


The results obtained using the above approximate formulation may be 
compared to the two previously given results and seen to be in satisfactory 
agreement. Graphical results, obtained through use of this last equation, 
were discussed at the beginning of this section. 
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3.0 SPATIAL FILTER EFFECT AND OPTIMAL FILTER FOR GEOID DATA PROCESSING 

3.1 Geoidal Data Processing Results 

The purpose of this section is to present derivations and results relating 
to optimum data processing methods and required data rates for CEOS-C altitude 
information. The first part of this section addresses the effects of (1) geoidal 
power spectral density, (2) spatial filter function, and (3) altimeter measure- 
ment random error on data processing characteristics. Weighting functions are 
derived for both undulation and slope filtering based on geoidal spectral data 
extracted from the Sky lab altimeter observations. The results are optimal in 
an additive Gaussian noise, minimum-mean- square error sense and largely represent 
an extension of the earlier work of Cohen and Zondek [1]. The second part of 
this section gives a detailed analysis of the spatial filter effect which results 
from the finite area illuminated by the altimeter. We assume that the altimeter 
spatial filter effect averages over all significant waveheight spectral com- 
ponents . 

Results of the analyses argue that the GEOS-C altimeter will be capable of 
measuring geoidal components to a short wavelength cutoff in the neighborhood 
of 20-40 km. for geographic regions containing pronounced short- wavelength 
features. For a satellite ground-track velocity of ~7 km/ sec, this corresponds 
to a Nyquist frequency of ~0.3 Hz. Taking this Nyquist rate and the time-expanse 
of the derived weighting functions as a measure of the geoidal information rate 
of the altimeter, the GEOS-C 10 sample/second data base is considered to be 
more than adequate. 

Figure 3.1 displays a power-spectral-density (PSD) plot for the Puerto 
Rican Trench region which was computed using fast Fourier Transform methods 
and a Hanning type convolutional window. The data base comprised SL-2 , Pass 4, 
Mode 5 with 100 and 130 nanosecond pulsewidths (pulse compression was not 
functioning during SL-2). The Puerto .R-± can Trench data was used since we desired 
to obtain PSD results for an anomalous region, which should contain more energy 
in short-wavelength components than anomaly-free regions. The PSD so obtained, 
and data processing results derived therefrom, should represent the best 
opportunity for the altimeter to obtain information relating to short wavelength 
undulations and should yield an approximate upper bound on data processing 
requirements. 



Power Spectral Density 



Fig. 3.1 Geo id Undulation Spectrum of Puerto Rican Trench Area and Wiener Filter 
Transfer Function. 
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Referring to Figure 3.1, the dashed line corresponds to the density level 
for which a 5 Hz rectangular bandwidth, white noise spectrum would yield an rms 
level equal to 0.5 meters. Since the S-193 tracker has approximately a 3.3 Hz 
equivalent noise bandwidth, this level consititutes an approximation to known 
GEOS-C intensive-mode system characteristics. The noise level shown in the 
calculated spectrum represents the Skylab altimeter noise level (1-2 meters rms) 
We will subsequently verify that the spatial filter function corresponds to 
considerably shorter wavelengths (less than 10 km) and that the calculated PSD 
is not contaminated by the altimeter footprint effect. 

Note that the observed spectrum represents an asymptotic behavior which in 

the frequency parameter (f) is approximately f . As discussed in [2], Kaula s 

-3 

jnodel of one- dimensional spectral behavior behaves as f 

Since observable geoidal components are of much longer wavelengths than 

spatial filter effects, we interpret the data in Figure 3.1 as a cut through a 

two-dimensional spectrum. In wave-number space (k ,k ) a directional spectrum 

x y 


S(k ,k ) with a k behavior will yield a one-dimensional asymptotic behavior 
of k (due to integration over the angular coordinate of the polar coordinate 
set). Therefore, we feel that the observed spectrum depicts the proper theo- 
retical behavior. 

The optimization technique we use is the Wiener-Hopf formulation, which 
for the correlation functions R(0 of signal s and observation y, gives the 
optimum' impulse response h Q (t) as the solution to the integral equation [3] 


R (T+r|) 
sy 


■/ 


h o (h)R y (T-y)dy x>0. 


For non-real time processing, an estimate of a value at time t can be based on 
both past and future values. Therefore, the proper lower limit on the integral 
is —CO and the integral equation becomes a convolution form which is readily 
solved by transform theory. ' For our purposes the form of the solution is [3] 


H(to) = 


S(w) 


S (w)+N(w) 


where S(u)) is the geoid undulation power spectrum and N(u>) is the additive noise 
spectrum. Since the altitude tracker has a noise equivalent bandwidth of ~5 Hz 
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and a random error standard deviation of "0.5 m, N(*) may be represented as a 

O'} 2 — 

white noise spectrum with a density of (.5) m/5 Hz = ,05 m /Hz or 7.96 x 10 

2 . 

m /radian. Using the break-point approximation (the asymptotes of which are 

shown in Figure 3.1) to S(u)) as 


S(aj) 


71.66 » 6.554 x 1Q~ 4 
w 4 -.0512w 2 + 6.554 x 10' 4 


the optimum transfer function is found to be 


V“> 2^ 

u -.0512W + 5.9006 

This function is also shown in Figure 3.1. Note that the asymptote is twice as 
steep as the spectral decay. At the intersection of the break-point spectral 
approximation and the GEOS-C noise level (which occurs at "22 km), H Q (w) 
introduces an attenuation of "12 dB. The 3 dB attenuation point occurs at ~40 
km. 

Figure 3.2 shows the spatial filter response function for the GEOS-C system 
which is derived in section 3.2. Note that the solution [H q (co)] given above 
effectively truncates geoidal data at considerably longer wavelengths than does 
the spatial filter effect (its 3 dB point occurs at "10 km). Had this not been 
the case, a considerably more involved Wiener-Hopf form would have been required 

The optimal filter H Q (u)) has been inverse Fourier transformed through use 
of contour integration, and the normalized impulse response found to be 

'h(t) = e“°‘ 8755t (cos 1. 289t + 0.6792 sin 1.289t) for t>0. 

j 

Knowing that the optimal geodetic slope filter is the derivative of the 
optimum undulation filter, the impulse response for slope estimation is 

-0.8755 e~°‘ 8755t (cos 1.289t + .6792 sin 1.289t) 

+ e“ 0#8755C (.8755 cos 1.289t - 1.289 sin 1.289t) for t>0. 

The undulation filter impulse response will be an even function of time, whereas 
the slope filter impulse response will be an odd function [1]. Both response 


d h(t) _ 
dt 
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functions are graphed in Figure 3.3 and tabulated values are given in Appendix 
A. This figure shows the overall characteristics of the two filters; for 
digital computer implementation the response functions must be truncated and 
properly re-normalized. To insure that the long wavelength structure is 
invariant under discrete-time filtering, the truncated convolution weighting 
coefficients should sum to unity. Also, in application of this form of 
filtering, the geoid power-spectral density of the actual area of interest 
should be investigated and the mean square error properties, given by 

CO 

. . Ij /»„<•>»<->-. 

.00 

should be compared with weighting functions which are computationally more 
efficient (e.g., rectangular or arithmetic moving average). It is interesting 
to note that the time expanse of an ideal rectangular impulse response function 
which has a sin x/x frequency domain response that matches H^Cu)) at the 3 dB 
point, corresponds to an averaging period or impulse response width of 2.22 
seconds. The h Q (t) form shown in Figure 3.3 has a width at the e ^ points of 
"2 seconds. 
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3.2 SPATIAL FILTER EFFECT 
3.2.1 Background Discussion 

The term spatial filtering effect as applied to radar altimetry denotes 
the ideal, i.e., noise-free, capability of the altimeter to resolve surface 
features. Spatial filtering is inherent in an altimeter system because of the 
non-zero width of the system point-target response. Thus, even if we consider 
the measurement process to be noise-free, the above factor limits the ability 
of the system to ’’map" small scale surface features. 

Typically, one calculates the radar footprint radius R for pulse-length 
limited geometries from the relationship 

R - y/hcT 

where h is the satellite altitude, c is the speed of light, and T is the radar 
pulse length. It is then assumed that the footprint diameter (2R) approximates 
the minimum surface wavelength which the altimeter can measure, and that the 
altimeter response is essentially flat for surface wavelengths greater than 
2R and zero for wavelengths less than 2R. In essence then, the altimeter is 
treated as an ideal (rectangular) low-pass filter whose cut-off frequency is 
determined by the transmitted pulse length. While such a characterization of 
the altimeter may be sufficient for "order-of-magnitude" calculations, it is 
inadequate for more precise system modeling. The problem of specifying altitude 
sampling rates and optimal data filtering require that the asymptotic (high 
frequency) behavior of the altimeter spatial filter characteristics be investi- 
gated. 

In a previous work [4], approximate formulas were obtained for determining 
the effect of a one-dimensional, sinusoidal, corrugated surface on the mean 

*w< m 

return waveform. Upon further analysis, we have determined that the approxi- 
mations used to obtain those formulas are not always valid. In this report we 
correct this formulation for the scattering process to obtain the effects of 
the surface on the mean return waveforms. For the GEOS-C altitude, pulse width 
and split-gate tracking- configuration we obtain results which illustrate the 
sensitivity of the system to surface undulation wavelengths. 

For the purposes of this report we will represent the geoidal surface 
undulation as a single sinusoid of small amplitude and very low frequency. In 
other words we will assume that on a very localized basis the geoidal perturba- 
tions of the mean flat sea may be represented as a sinusoid of relatively small 
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amplitude and low frequency. The analysis could be extended to the case of two 
or more surface harmonics; however, the analysis and computations become 
prohibitively tedious. 

The geometry of the problem is shown in Figure 3.4 along with a definition 
of pertinent terms. We assume that the surface undulations vary in one 
direction only. That is, the deterministic component of the surface height may 
be represented as 

z = a sin(k pcos4> - ot ) 
s s s s 

where a is the peak amplitude of the undulation, k g = 2TT/ is the wave-number, 
and a is the "phase" of the surface directly beneath the altimeter. For the 
time being, the altimeter will be considered to be fixed relative to the surface. 


3.2.2 Determination of the Mean Return Waveform 

The mean return waveform is given by the following expression [5] 


P (t) = 
r 


(4tt) 


/ 


P (t- %)G 2 (Q,u>)g 0 (Q,w) 
t u _ 


dA 


( 1 ) 


SCATTERING 

AREA 


where t is the time relative to the time of transmission and dA is the elemental 
scattering area on the ocean surface. For purposes of this analysis, we consider 
the altimeter antenna to be nadir pointed, thus w - 4> and 0 = We also will 
only consider the case of pulse width limited geometries and very short pulses. 


i.e. 


P (t) s 

r 


/ y t- % 

(47 


( 2 ) 


where G is the boresight gain of the antenna and o> is the surface scattering 
o 

cross-section per unit scattering surface at 0 angle of incidence. If we 

define r as the distance from the altimeter to the point on the z = o surface 
o 

which is at the same radial distance p as the intersection point of r with the 
true mean surface, then from Figure 3.5 


r 


r 

o 


z 

s 


V l+Cp/h ) 2 





ELEMENTAL SCATTERING 
AREA 



SOMEWHAT LARGER THAN THE RADAR "FOOTPRINT," 
THE SURFACE WAVE NUMBER k IS k s IN MAGNITUDE, 

IS IN THE x-DIRECTlON, AND HAS A PHASE ANGLE a* 
RELATIVE TO x ■ 0. THE AMPLITUDE OF k s IS a„ SO 
THE SURFACE ELEVATION AT ANY POINT IS GIVEN 
BY z s , 

z t = a s sin(k s pcos0 - a,) 


THE SUBSATELUTE OR NADIR POINT IS THE 
ORIGIN OF THE x,y,z OR p,0.z COORDINATE 
SYSTEMS. 

$ = ANTENNA POINTING ANGLE, OFF NADIR. 

0 * ANGULAR LOCATION OF SCATTERING 
AREA RELATIVE TO BORESIGHT AXIS. 

5 * POLAR ANGULAR LOCATION OF BORE- 
SIGHT POINT RELATIVE TO x-AXIS. 

0 » POLAR ANGULAR LOCATION OF 

SCATTERING AREA RELATIVE TO x-AXIS. 

p - POLAR RADIAL LOCATION OF 
SCATTERING AREA. 

= PROJECTION ON z = 0 PLANE OF THE 
ANTENNA PATTERN'S AZIMUTHAL 
ANGLE. 

k, * 2 jt/X s . 

h “ SATELLITE ALTITUDE. 

r - DISTANCE FROM SCATTERING AREA TO 
SATELLITE. 

r « hKl +(p/h)2-22,/h 


Fig. 3.4 Satellite-Ocean Surface Geometry and Summary of Notation 
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where r = Vp +h and z = a sin(k pcos4>-cx ); thus 


s s s 


t 2r _ 1 

c c c Vi+(p/h? 


where c is the speed of light. For ^«h. 


_ nt fJl 2z 

t - ^ ^ - % + -*-■ 

c c ch c 


Equation (2) becomes 

rJZ 2z 

A (t .as. 

—4 — 1\ , , _ 2h 

where we have approximated the r factor bt r Q . Substituting T - t — q 

and approximating dA by pdpdcj), i.e., we integrate over the z = 0 plane but 

retain the z g dependence in P^, yields 


-£.■+ 


A 2 G 2 a° r Y t (x ” he + c 5 

p (t + — ) = — 2=- / t ?apa<t>. 

r C (4ir) 3 J 


We now take the transmitted pulse to be Gaussian, i.e.. 


P t (t) = P T e 


then (4) becomes 


r c n (4 tt) j J J 


„ , ~2 2z 2 

” 2lT - -i- (t - fi- + — 

f [ e - 2 hc c . 


(l+^/h 2 ) 


2 

Substituting (S' - in (5) yields 


2 2 ' 

p f + 2 h. _ XGgP T 
P ( T+ ) ~ 3 3 

r c 2ti (4 ti) 


00 2tt 


- J ^exp | — ~(p-T“2z s /c) 2 | dpd^ 
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~2 

where we have ignored the factor (1+cp/h) since it has little effect relative 
to the Gaussian dependence on p. With the change of variables in P, z g becomes 


z - a sin 
s s 



After some trigonometric manipulations, the exponent inside the integral in 
(6) becomes 


4a 


n o i» " 

[p-T-2z /e] = (p-t) - (p-T)sin(k yhcp cost}) - a ) 



[1 -cos (2k \hcQ cost}) - 2a )] 
S 1 s 


(7) 


. We then use the following identities [6] to remove the last two factors from 
the exponent: 

2a s 

— £-(p-T)sin($cos<Ji - a g ) 
ccr 


= I 0 (x) 



+ 2 (-l) m I 2 m +i^ x ^ sin|(2m+l) (gcos<|)-a s ) | 


m=o 


+ 2 'y ^ (-l) m I 2m (x) cos|2m(gcoscJ)-a x )| 


m-1 


( 8 ) 


and 


exp < - 


2 2 . f 

a a ) a e 2 ( 

S + ~ 7 ~ cos ( 23 cos<J) - 2a D )[ = exp -(— ) <I Q (y) 


2 2 2 2 
c a c a 


f 


where 


+ 2 I n (y) cos |2n(£ cos<}> -a g )jj 


n=l 


3 = k \/hcp 


(9) 


2a 


ca 


-<P-t) y 


2 2 
c a 


x 
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and the I (•) are Bessel functions of the second kind and order n. Taking the 
n 

product of (8) and (9), expanding the trigonometric in series of Bessel 
functions of the first kind and integrating terra by term yields the following: 


2tt 


( J 1 

/ exp j- — 

J l 20 


2 z 2 ) 

2 [p-T- -£-] [ d<|> = 2 tt e 


( V 


20 


2 v c a J 


i 0 (x) x 0 (y) 


+ 21 (x) 
o 


2 r n (y>J o 


(2n$)cos (2ncx g ) 


n=l 


- EE (-l) m * 2n ri-l W X n (y)6 n o [ J o ( [ 2m+1 - 2n l 3) si n ( [ 2nri-l-2n ] a fl ) 


m=o n=o 


+ J o ( [ 2m+l+2 n ] 3 ) s in ( [ 2m+l+2n ] ct g ) J 

-EE (-1) ^2m^ x ^n^^n,o 3)c° s ( [- ^ 

m=l n-o 

+ J ([2m+2n])3)cos([2m+2n]a ) 
o ° 



where 


6 . 
n,o 


1/2 n = o 
1 n = 1,2, *• • 


Rearranging and changing the indices on the various terms in the above leads 
to the following result; 



• N ' ^ I (x)I (y)S 6 ) J ([m-2n30) 

/ j / ^ m n J ' n,o m,o o 
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^l+(-l) m }cos^cos([m-2n3a s ) - j l-(-l) m [cos (m-l)^sin( [m-2n]a s )J 
+ J o ([m+2n]3) ^Jl+(“l) m }cos^cos ([m+2n]ct s ) 


[l-(“l) m ]cos[ (m-l)~] sin([m+2a]a ) 

s 


( 10 ) 


Equation (10) is the result of integrating (6) with respect to the ^variable; 
it is exact in that no approximations were employed to obtain (10) from (6). 
The right hand side of (10) may be simplified considerably when 

2 

* = 77 ~ 

or the pulse length is much larger than the peak amplitude of the sinusoidal 
undulation. Under this condition we note that 


I n (y) ~ T 1 * 

2' r(n+l) 

and we only need to take the first two terms in the I series. Thus, the mean 

n 

return waveform is given by 


9 , A 2 G 2 g°P ire - 
P (x+ Ih) s _o T 

r c / . % 3, 3 

(4tt) h 


a_ 2 1 . 2 

“7 [P“ T 3 
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J e 


20 


o o 

xl (^ P " Tl ) a B,o | [ (_1)m+1 3 c0s ^[ J o ( n ^ s v / ^cp)cos(ma s ) 

m=o u 


+ 1 J o ^ [ m ' 2 3 k s Vhcp)cos ( [m-2 ]a ) 


+ J 0 ([n^2]k s Vh^)cos([m+2]a s )[] -[l-(-i) m ] CO s^l)3L 


[ J 0 ( mk s VTicp)sin ma g + |(“) 2 |J Q ([m-2 jk s V^)sin([m-2]cx s ) 
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+ J ([m+2]k Vhcp)sin([m+2 ]« )i >dp< 

O S SI | . 


( 11 ) 


Equation (11) represents the mean return waveform scattered from a unidirectional 
sinusoidal surface under the following assumptions: 

(a) pulsewidth limited geometry, and 

(b) the peak amplitude of the surface undulation is small compared 

a 

to the pulsewidth, i.e., << O'. 

Expanding the first few terms of (11) yields 

a 


.2 2 OTi “ s .2 /* 

_ , . 2h, X G o° V c - ( W> / 20 

p r (x + — ) rr e J e 

r c (4 tt) ri J 


2 [ T ] 2a 


I 0 (-^[P-T]) 

ccr 


* [l + (cf) 2 V 2k sV^ )cos 2a 1 -2 I 1 (— |[p-T]) 

L J CQ 

[ Jo (k s V^P) S ma s + ^(~) 2 |-J 0 ( k s Vhcp) sinCi s + J Q (3k s Vh^p)sin 3a s |] 

2a r* -t *2- o ( 

-2X 2 (— |[p-x]) J o (2k sA /hcp)cos 2a g + 1 

ca L ° ( 

+ J q (4k g hep) cos 4ajj + 2I 3 ( — ~[p-T ] ) |j Q (3k g yj hep) sin 3a g 


* a s + J Q (5kg "^icpjsin 5a & | 


+■** dp. 


Of all the terms appearing in (12), only the first can be approximately 
integrated in closed form. For fs. «T*a, as shown in Appendix B., 
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00 1 2 
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^ *^2rra 


(-S.) 2 a 2 - 
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( 12 ) 


Thus, 


p(T + 2h )= ^Vi 

p r^ T + c } ,, N 3, 3 
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C > + VTT> 
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a 2 00 
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1 <-*• “ r 

J Q (3k g v1icp)sin 3a g + ^C") d o (k g v4icp)sina 


4'J o (5k g /hcp)sin5aj | 


and since 


dp 


a 2 o 

-(— ) a 2 

> ca erf (— — -) + I (-S 5 ,) ~ erf (-*-—) + 1, 

o V 2 C O a Si 


(13) 


(14) 


the first square bracketed term in (13) is essentially the flat-sea mean 
scattered return. The term in (13) involving I 2 (^f[P- T l> which is independent 
of k is a higher order correction to (14) which may be ignored when a s /c«a. 
The remaining terms in (13) depend on k g and thus exhibit the dependence of 
the mean return waveform on the surface wave number. Provided k g v^ic* is not 
too large, the dominant term inside th'e' integral in (13) is 


-2sina. 


/ 


- -^[P-T] 2 
20 


2a 

I l < ‘ — | [p “ T ^ J o^ k s »^P) dp * 
ca 


The exact value of k She for which the above term no longer dominates the k s 

s 

dependent terms in (13) is not analytically obtainable, but may be readily 
determined by numerical integration of (13) . 
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For a flat sea, the mean return, waveform is a convolution of the five 
following factors: 

(a) transmitted waveform, 

(b) flat sea impulse response, 

(c) radar observed waveheight distribution 

(d) radar receiver impulse response, and 

(e) tracker loop time positioning jitter. 

However, for the non-flat sea surface such as we are dealing with in this 
example, it is not possible to deal with an "undulating sea impulse response" 
or the analog of (b). That is, if we found the response of the undulating 
surface to an impulse function, the response to another pulse shape would not 
be equal to a simple convolution. This can be seen by referring to equation 
(1). The convolution property arises because for a mean flat sea, dA=rdrd$ 
and the <J>— integration can be accomplished independent of r and the argument 
of P t remains as the difference t-2r/c. For the undulating sea dA^rdrd^ and 
r is a function of both <f> and (S'; thus, after performing the (^-integration i n 
(1) the argument of P t will not necessarily depend on the difference t-2r/c. 

For this reason, the concept of a "undulating sea impulse response" has no 
meaning. In fact, the concept of linear scatter theory or the entire multiple 
convolution model of the process may be questionable; for conditions under which 
geoidal and ocean surface wave lengths are of comparable length. However, this 
topic is much beyond the scope of this investigation. 

For the purposes of this report, we will assume the following: the system 

point target response* of the altimeter is Gaussian with a 3dB pulsewidth of 

T ; the radar observed waveheight probability density function is Gaussian with 
P 

an rms waveheight equal to 0 gg and the tracker loop time positioning jitter is a 
continuous Gaussian process with standard deviation equal to Cj . If we assume 
that the significant waveheight of the waves (H^ ) is equal to four times the 
rms waveheight, the pulse width parameter 0 is given by 


V<- 


508 H 3/ 3 )2 + ( * 602 T p )2 + a j 


where H-k has units of feet and T and a. have units of nanoseconds. For the 
^3 P 3 

GEOS-C altimeter system * 3 ns. For this study we have chosen a composite 


*The point target response of the system is the convolution of (a) and (d) , above. 
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a of 8.16 ns. Table I shows the combinations of T and Hi which satisfy the 

p V 3 

above equation with o = 8.16 ns and *= 3 ns. If the system point target 
response pulsewidth is 12.5 ns then Hi, = 1.98 feet which is a small sea. 

3.2.3 Determination of the Split-Gate Tracker Response 

Equation (13) represents the mean return waveform as influenced by a 

random scattering surface having a mean sinusoidal profile. In order to 

determine how the altimeter tracker responds to this return waveform, we must 

apply the split-gate tracking algorithm to the return. We first normalize 
2ti 

P r (l + — ) such as would be done by the receiver AGC, i.e.. 


P r (T) = 


P (T + 2h/c) 


where 


r = 


,2 2 

X G P^TTC 
o T 

(4it) 3 h 3 


We now integrate over the ramp part* of the normalized return to form the 

function H„(k ,a ,a ) where 

R s s- s t 2 


WW = J 


P r (T)dT, 


(15) 


-T /2 
P 


TABLE I 

Combinations of pulsewidth and significant waveheight which yield o = 8.16 ns 
with = 3. ns. 


T 

P ' 

Pulsewidth (ns) 

h i/ 3 

Significant Waveheight (feet) 

8 

11.55 

10 

9.1 

12.5 

» 

1.98 


*For GEOS-C the tracking gates have the same width as the system point target 
response, i.e., T p . 
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while the plateau gate function H (k ,a ,a ) is formed by integrating P ( T ) 

p S S S i 

over the same time interval but starting at a time well into the plateau 


region of the mean return, i.e., 


T+T 


H (k 
P s 


,a ,a )= / 

a s J 


_p/2 

P (T)dT, 


(16) 


T-T 


P/2 


where T » T . 

P 

When a = o, the normalized mean return becomes 
s 


lim P ( T ) = [1+erf (— ) ] , 

a s^° r 2 a /2 

thus 

T 

\^s> a s = °* a s > = ^ cr(-^) (17) 

and 

H (k ,a = o,a ) = JZtF a T (18) 

p s s s P 

Comparing (17) and (18) we see that 

V k s ,a s “ “’V = 2 H p°V a s = °> a s ) (19) 

and, thus, for the mean flat sea, the tracking law is satisfied when the 
integrated ramp and pleateau gate values are defined as in (15) and (16). 

In the actual system the altitude tracker develops an error signal for 
each received pulse, each of which is an ensemble member of the mean power 
versus time- relationship of the backscattered signal. Averaging characteristics 
of the closed-loop system thus provide. a measure of mean signal properties. 

Since we can only -calculate signal statistics - not individual ensemble members 
it is necessary to conceptually reverse the order of the altitude measurement 
and averaging processes. This amounts to an ergodic assumption, which previous 
Monte Carlo simulation studies have shown to apply. Also, since the altimeter 
uses a square- law detector, the averaged video signal should correspond to the 
calculated P r (t) for the noise-free case. 

We further assume that the mean return waveform as a function of k does 

b 

not differ appreciably from a shifted replica of the mean flat sea return and. 
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also, that 2a /c « a, i.e., the maximum possible two-way time shift due to 
the sinusoidal mean surface is small compared to a pulse length. Under these 
assumptions we can mathematically replace the closed-loop tracking process by 
a fixed discriminator relationship, i.e., 

e t = ejy'vvy - 1 y vw ] <20) 

where e ^ is the time position error relative to 2h/c and we define 3 such 

that when k = 0 (X = °°) and oc — tt/ 2 or the surface is flat but a meters 
s s s . a 

below the e. = 0 sea surface then e = 2a /c. 

t r s 

To obtain (20) we note first of all that if T is sufficiently large (as 
defined in equation (16)) then 

H p (k s’ a s ,a s ) * y k s’ a s =O,0t s ) = S2* aT v > 

or changes in the location of the leading edge of the pulse as a function of 

k have no effect on the plateau gate integrated value. Thus, (20) reduces to 
s 

VWV - ~r 0 V • 

When k = o and a = tt/2, we determine 3 so ^that e. = 2a /c, i.e., 
s s t s 




2a s /c 


/2tT 


Hr. ( k s = ° * a s ,ot s =1T/2 } r° T 


J 


When k = o, the normalized mean return waveform may readily be obtained from 
s 

equation (6) and 


, . - / . ' 2h, /2ia 

lim P (t + -r) - — o — 
, r c ^ 

k -*o 
s 

2h 


/T + 

1 + erf( 


^ a s sin (-a ) 
c s 

a J2 




Integrating lim P (t + — ) between the limits of -T 12 and +T /2 and setting 
° r , c P P 


k s +o 

a = tt/ 2, we find 
>s 


H R (k s =o,a s ,a s = 1 r/2) = {T p +(~f - -^)erf 


D 



47 


T 2a - T za rt /T Orr 

, p . , 1 r P . s 1x , O y2 20 

- (n + )erf ( [Hr + J ) + e 

2 c 


T 2a 

i + — 

2 c 


, T 2a 2 
_ J_ ( _E _ _JL > 


2 2 c 


a/2 ‘ /? 

, T 2a 2 

a ^ 2a 2 2 c 

■ e 

/ir 


> - 


( 21 ) 


Substituting a = 8.16 ns, T = 12.5 ns and a = 0.4 meter in the above yields 

P s 

vv°-w* /2) - J iF L T P = - 29 - 898 > 

and 

3 = -29.898. 

The time shift of the tracking loop is therefore given by 


e 


t 


(2a s /c ) f 
29.898 



( 22 ) 


where H^Ck ,a ,ct ) is 8 iven by (15) with P r (T) computed from (13). 

The height of a point on the surface directly beneath the altimeter is 
given by 

h A = h “ Z s 


or 


h. = h-a sin(k pcos<4-a ) 
ASS S 


(23) 


where h is the height of the altimeter above the mean flat sea. The phase ot g 
can be put into a one to one correspondence with the spacecraft velocity in the 
x-direction by the following: 


where V x is the spacecraft velocity in the x-direction and t is spacecraft 
time. Setting £=o in (23), results in the following altitude profile encountered 
by, the altimeter as it moves with constant velocity in the x-direction: 

* h. = h + a sina (24) 

Ass 

or 

h. - h + a sin(k V t). 

A s s x 
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From (24), we see that the time displacement, e t , of the tracker should be 
proportional to the sine of the surface phase angle. As noted previously, for 
relatively small values of k s /he such will be the case; however, when k g /he" 
becomes large, the dependence of on a g is no longer sinusoidal. 

To demonstrate the sinusoidal behavior, we return to equation (13) and 
compute the AGC’d mean return in the limit of k^+o. The details of the 
integrations are given In Appendix B and 

2 

/2?c t , . a , _/?¥ Jj > — o„ 


Xim P (T + % 

V° 


|erf(— - — •) + lW 

l rr VT ’ 


a / 2 


cos 2a_{l + erf ( 

1 co 


T —)\ 
a A ' 


- * (-f> 


T 2 

t 1 " a/I }][ sina s + f si 


+ sin 3a 


}]" 2 l 2 a t + a 


1 + erf(— - — ) - 


4t 


/2a 


a SZ /tt 


>] 


. a 2 

1 , s 


• fcos 2a + tt(— ~) (1+cos 4a )] +• • * 
l s z ca s 


and T 

J2. 
2 


/ 


-T 

_E 

2 


lim P (t + — )dT = ° (T ) + 

k **o r c 2 p 

s 


/B a 


2 r cos 2 Vy 

c a 


- 2 ef)(T p )[slna s +.|(^§) {-sina s + sin 3a s }j 


a 2 SB r 1 a 2 

2 — ( T p ) cos 2a s + 2^co^ ^ + cos 4a s ) I + ’“* 


2c a 

therefore the tracker 'loop time positioning becomes 

2 


-] 


, lim G 
k -K) 
s 


t - - T p { ■ Jr cos 2 “s - 2(! f> [ si ™ : 
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a 2 

SI 

+ ) < -sin a + sin 3a, 


a 2 fit 


2W {' sil 


a v iii r 

/ J 2c a L 


~ ,cos 2a 

4. . I S 


+ |(J) (1 + COS 4a s 


>] ...| 


Substituting o - 8.16 ns, a g = 0.4 meter yields 


<2a s /c) 

k'+o E t = 29.898 

s 


T <.273 cos 2a - 2.67[sina + (.013)(-sina 

p ( S'. S s 


+ sin 3a ) ] - .273[cos 2a + (.013) (1 + cos 4a ) ] + 

s s s ^ 


or 


2a 

,lim e = (1.11) (— %sina . 

ic -K* t C S 

s 


(25) 


(26) 


The factor (1.11) is due to the fact that the integration of (13) as shown in 

Appendix B is approximate while the denominator in (25), i.e., 29.898, was 

exact (see equation (21)). The important point to note in (26) is that the 

tracker loop time positioning is proportional to sin a g . 

When k is not equal to zero, it is not an easy task to demonstrate the 
3 

dependence of e on sin a g . Because of the complexity of the integrand in 

(13), it is not possible to analytically obtain the value of k g for which the 

terras multiplying sina are dominant. That is, we would like to find the value 

of k such that for k < k 
s s s 


“ , . 2h N a 

P r (T + T 5 * 2~ 


a 2 0 

-(— ) a 2 ‘ 

> 00 erf(— *-) + X Jrri> 

a / 2 ° c 2 a 2 


a 2 Cp-T) 

-(— ) r 2a T i a o 2 1 ?rr 2 

- 2 sina_ e C ° / I (— §[p-T])J (k /hcp)|l - j(“) Je 


s 


dp 


(27) 


and for e we then have T 

a 2 ~ 

e. » + 1 2 sina g e'^ f * J £f[ P-T ^ (k, /h^) 

-T o 

2 


(2a /c) 


29.898 


c o 
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• [■ - K§> 2 ]*’ “ 2> •»«} • 

Equation (28) can be rewritten in the following form: 

E t “ (-^)A(k s )sina g 

where A(k ) has a maximum value of one when k — o and it determines the 
s s 

dependence of £,_ on k . Because of the J function in (28), A(k ) will go to 
t s o s 

zero for certain values of k and as k increases beyond this value, A(k ) 
will change sign. A change in the sign of A(k g ) indicates that while the 
altimeter may be profiling the gross variation of the surface, i.e., sin ct g , 
properly, it is interpreting peaks and valleys in the surface as valleys and 
peaks, respectively. With the limited amount of computation that we have 
accomplished on this problem, we have found that (28) is a reasonable approxi- 
mation to (22) for k g slightly less than that value for which A(k g ) goes to 
zero for the second time. For values of k g beyond this point, the variation of 
£ with surface phase is no longer sinusoidal since other ot dependent terms in 
(13) become dominant. However, for spatial wave lengths at which (28) becomes 
invalid, an absurd noise level for altimeter observations would be required 
to permit compensation of spatial filter effects. 

3.2.4 Results for GEOS-C Intensive Mode of Operation 

Much of the preceding material has been concerned with the spatial filter 
problem in general. In particular we have attempted to show how the sinusoidal 
surface affects the mean return waveform and how this effect may be translated 
into the profiling capability of the altimeter’s tracking loop. We have also 
pointed out some of the simplifying assumptions that can be made in dealing 
with the rather complicated form of the mean return waveform. We will now 
address the GEOS-C problem. As before we will assume a 3 dB pulse width (T ) 
of 12.5 ns for the system point target response, a significant waveheight of 
1.98 feet and a tracking loop jitter standard deviation of 3 ns. The orbiting 
altitude of the altimeter will be taken as 880 km, while both the ramp and 
plateau tracking gates will be assumed to be 12.5 ns wide. The system parameters 
assumed represent our present knowledge of the expected final configuration of 
GEOS-C in the Intensive Mode while the altitude and waveheight parameters are 
taken to be nominal operating conditions. 


( 28 ) 
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For the parameters listed above, we have computed the mean return waveform 
by numerically integrating equation (12) including all terms out to I^(2a s [P-x]/ 
ca 2 ). All computations were for a peak surface amplitude of 0.4 m. The 
resulting waveforms were integrated over the time interval - '-E- < T < -2. and 
the output of the tracking loop time discriminator was computed using (22) . 

A plot of e t , normalized by 2a s /c and as a function of surface wave number and 
wavelength is shown in Figure 3.6. This plot illustrates the extent to which 
the profiling ability of the altimeter degrades as the wavelength of the surface 
undulation decreases. Due to sea-state and jitter effects, the equivalent 3 dB 
pulse width of the system increases from 12.5 ns to 13.55 ns. The approximate 
cutoff wavelength of such a system is given by 


X c = 2-^/hc (13.55) 

or X =3.78 km. As shown in Figure 3.6, A is slightly smaller than the first 
c c 

zero in e . As A decreases beyond the point for which £ = o, £ changes sign, 

t s L 

That is, the altimeter indicates the presence of a valley in the surface when 
it is actually measuring a peak. Near the second zero of the altimeter 

profiling capability entirely breaks down because the altimeter no longer sees 
a sinusoidal surface. Thus, for A g < 2 km the altimeter is no longer profiling 
the surface. Such an observation can not be made by just examining Figure 3.6; 
this statement requires that the variation of £ fc with surface phase, a g , be 
examined and compared with the sina s< Thus, the plot in figure 3.6 should be 
"cutoff" at X g - 2 km. 

While Figure 3.6 is the desired output of this study, it is also interesting 
to examine how the mean return varies as a function of surface wavelength. The 
solid line in Figure 3.7 is that portion of the normalized return for a mean 
flat sea surface, i.e., = o. Thedashed curve for X^ = 31.4 km is essentially 

a shifted replica of the flat sea curve. It is shifted to the right (later in 
time) because the nadir point on the surface is a "valley" in the sinusoidal 
mean surface, i.e., a sin(-ct ) = -a . As the wavelength of the undulating 
surface decreases, the mean return shifts closer to the flat sea curve and also 
begins to change shape’. The curve in Figure 3.7 for X g = 6.98 km corresponds 
to the -5.5 dB point on the £ t versus X s curve in Figure 3. Figure 3.8 illustrates 
how returns from surfaces having — 4.49 km and 2.73 km compare with the flat 
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Fig. 3.6 A Plot of the Altimeter Tracker Response to a Sinusoidal Surface as the 
Wavelength of the Surface Varies, 
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sea return. For the 4.49 km return, the portion nearest the plateau region 
actually appears earlier in time than the flat sea. For k^ = 2.73 km, the 
entire portion of the return shown in Figure 3.8 appears earlier in time than 
the flat sea return, thus giving the appearance of having been scattered from 
an elevated surface. This character of the return, of course, gives rise to 
the sign reversal in £ t , i.e., the altimeter interprets peaks in the surface 
as valleys. 
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Appendix A 


This Appendix contains a tabulation of weighting coefficients for the 
optimum undulation and slope extraction filters. The n T n column is 
incremented by 0.1 seconds corresponding to the output altitude data rate. 
The "H" column is a list of the corresponding weights for extracting geoidal 
undulation data. The. "DHDT" column is a list of the weights for extracting 
geoidal slope data. As discussed in section 3.1, both columns of weights 
should be renormalized based on the truncation point selected for the data 
processing. 
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Appendix B 

In the course of the analysis presented in section 3.2, it will be 
necessary to evaluate integrals of the form 


■/ 


- 2a- 


-<P-t) 


I »(^ [p_Tl ) dp 


(Bl) 


I = I e 
o 

The purpose of this Appendix is to demonstrate how these integrals may be 

evaluated. First, we substitute 7] = p-T in (Bl) , then 

2a 


s 

/ - 2 o £ /2a \ 
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dn 


-T 


which can be split into two integrals, i.e. 


/ 2.0 /2 a v r 

6 I n(^ n J dn + J e 1L l — f r\ J dn 

0 -T 

The first integral can be found in standard references on integrals of 
Bessel functions and 
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n ( 11 ) 
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For the second integral in (B2), we will only be concerned with that portion 

of the return for which I T | <_ T . *, thus 

s/ « 


2a 2a /T 

-f hi < -f S 

ca z ca z V 


4 )- 


Because of jitter and sea-state effects, the 3dB pulsewidth of the system 


*T g is defined as the 3 dB width of the system point target response. 
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point target response will always be less than or equal to 1.67a, thus 
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ca — ca ca 


For the case to be considered here, a - 0.4m and a = 8.16ns, thus 
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Therefore, for the second integral in (£2) we may use the small argument 
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The above integral can be evaluated by noting that 
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If we also use the fact that( S / _) 2 << 1 , then 
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and substituting the above in (B3) end using (B.4) , (B5) and (E6) » we have 



or after combining terms 



0 
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If we substitute the values of o (8.16ns) and a (0.4m.) in the above 
coefficients of the bracketed terms we find that 

^ g = 10.22 

a 

— = 1.33 

c 

a 2 /27i 

— « 0.136 

4c 2 a 

Thus, each coefficient decreases by an order of magnitude . from the prior 
value . 



64 


4.0 ALT ITUDE- WAVEF OFM BIAS EFFECTS AND POINTING; ANGLE ESTIMATION 

4.1 Waveform Bias Effects on the Altitude 

The GEOS-C radar altimeter employs a split-gate time discriminator to 

extract altitude error on a pulse-by-pulse basis. A typical split-gate 

arrangement is shown in Figure 4.1 along with an idealized "flat-sea" average 

return waveform. A discriminator curve is generated by integrating the return 

waveform over the time extent of the ramp and plateau gates, multiplying the 

plateau gate integrated value by 0.5, and subtracting this value from the 

ramp gate integrated value. A typical discriminator curve is shown in 

Figure 4.2. The discriminator curve is a function of the time variable t^ 

which is defined (in Figure 4,1) as the shift in position of the two gates 

relative to the mean return waveform. For purposes of this study, we will set 

t = 0 when the center of the ramp gate occurs at a round trip delay time of 
g 

2h/c where h is the true altitude of the altimeter above mean sea level and c 
is the speed of light. Since the ramp and plat-eau gates are separated by a 
fixed time interval T, a shift in the ramp gate also results in an equal shift 
in position of the plateau gate. It should also be noted that when the center 
of the ramp gate occurs later (or earlier) in time than 2h/c, t will be 

O 

greater (or less) than zero. 

The three most important characteristics of the discriminator curve are; 

(1) the width of the linear position of the curve, (2) the slope of the curve 

about e(t ) =0, and (3) the shift (or bias) of the intercept or e(t ) = 0 
S'. g 

from the t =0 point. A large linear range of the discriminator curve 
S 

is desirable since it permits a one-to-one correspondence between e(t ) and 

O 

t for large tracker excursions. The slope must be selected based on sensitivity 
6 

and stability criteria which are dictated by the design of the remaining portion 
of the tracker loop. The bias or shift of the e(t ) = 0 point away from t = 0 
is a consequence of changes in the mean return waveform [1]. All of the above 
three characteristics are functions of the mean return waveform. In this 
study, we will investigate how the linear range, slope and bias of the GEOS-C 
discriminator curve depend on pointing error and sea-state. We will also 
investigate a slightly different gate configuration which results in a reduced 
bias for small or moderate waveheights without any sacrifice in the discriminator 
curve’s linear range or slope. This alternate gate configuration demonstrates 
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Fig. 4.2 Typical Disciminator Curve Resulting from a Split-Gate Configuration. 





67 


the desirability of an adaptive tracker system. No consideration was given 
to a 3 or 1 ns system since the results presented here for a GEOS-C type system 
indicate that the bias errors will be very small and, in addition, knowledge 
of the "sea-state” will permit correction of the bias to a very small value. 

4.1.1 Analysis 


The mean return waveform for short-pulse altimeters operated at or near 
nadir is a convolution [2] of the system point target response, the flat-sea 
impulse response, the waveheight probability density function, and the tracking 
loop jitter. For purposes of this study, we will ignore the tracker loop time 
jitter since it will be relatively small for GEOS-C, i,e., " 3 ns. A 

Gaussian function will be used to describe the system point target response 
■since previous studies [3] have shown that this is a reasonable approximation. 
That is, if we denote P as the point target response, then 


P pt (T) - e 


“ <2 V 


CD 


and the half-power pulsewidth (PW) is related to the standard deviation of the 
Gaussian approximation (Op) by the following, 

a 2 = 0.181(PW) 2 . 

r 

Similarly, the waveheight probability density function will be assumed to be 
Gaussian, i.e., 



P ss<*> 




( 2 ) 


where 0 is the rms roughness of the ocean surface. The convolution of (1) 
and (2), with z converted to two-way ranging time, yields the following 
Gaussian function 


P (t)*P (z) 

pt ss 
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n 2 2 
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where 


a 



2 

■) + 1 . 


For the Intensive Mode of GEOS-C, we are dealing with a pulsewidth 
limited geometry, thus we need only consider the shape of the antenna pattern 
very near bore sight in computing the mean return waveform. We may, therefore, 
approximate the antenna pattern as circularly symmetric and Gaussian, i.e.. 


2 • 2 fl 

G a ( 9) - G o e - Y® ln ® 


(4) 


where G is the boresight gain, 0 is the angle measured with respect to 
o 

boresight, and Y is the antenna pattern taper factor. If we match (4) to the 
measured pattern at the 1/2-power point, we find that 


Y “ 2.895 sin 2 (^|), 


(5) 


where BW is the 3 dB beamwidth of the measured antenna pattern. Under the 
Gaussian antenna pattern assumption, it has been previously shown [3] that the 
flat-sea impulse response is given by 
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(6) 


for T > 0 and P^T,^) “ 0 for T < 0, where £ is the pointing angle of the 
antenna boresight with repsect to nadir. The complete mean return waveform 
is convolution of (6) and (3), i.e.. 
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and the I (*) are Bessel functions of the second kind and order n. If we let 

A ^ 

— *T j Ct At v /2 and q = ca ct/h /2 then (7) becomes 
P P 
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( 8 ) 


For GEOS-C, the point target response half-power pulsewidth (PW) will -be 

assumed to be 12.5 ns;. .thus, a = 5.32 ns. Since the beamwidth of the GEOS-C 

P -3 

antenna is about 2.6 degrees, Y = 1.49 x 10 . The altitude of the GEOS-C 

spacecraft will be about 880 km, hence q = 1.28 x 10 6 . Provided we restrict 
the range of T to less than, say 500 ns, we can ignore all terms in the 
series in (8) except for 1 Q (*) since for £ £ 3° the factor raised to the n — 
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power is very small but dominant. To evaluate the remaining integral, we note 
that all terms inside the integral show a small variation compared to the 
Gaussian pulse shape factor, and also 

1- (qx+1) 2 cos^£ ~ l-cos 2 £(l-2qx) 

2 2 

= sin £ + 2qx cos £ 

and 

[1- (qx+1) ^]sin 2 £ ~ 2qx sin 2 £ 


/(qx+1) -1 
(qx+1) 

Hence, for 0 < x £ 500 ns and £ < 3° , 

( 4 2 2a I 
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When x < 0, we see from (8) that the correct expression for P r (x,£) is 
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The approximate expressions given by (10) and (11) have been checked with 
results obtained by numerically integrating (8) and the agreement was found to 
be excellent. 

We normalize (10) and (11) to account for the action of the AGC and form 
P r (T,0 where 


P„(t, 0 = A(x) |i[l+erf( ^ )]} , 

r ' a a *2 ' 

P 


and 


4c 


, . Y h cos2 ^ T , 4 Yc ^ 0 


A(r) - 


I (- 

° Y /h 


T < 0. 


( 12 ) 


(13) 


The discriminator curve is obtained by integrating (12) over the time expanse 
of the ramp and plateau gates, multiplying the integrated plateau value by 
0.5 and subtracting this value from the integrated ramp value. If E (tg) is 
the discriminator curve, then 


T +t s 


T + 
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-T 


f +t g 
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+ t 


where we have assumed that the ramp and plateau gate lengths are equal and 

T g = gate length, 

w(T) = gate weighting function, 

T = time delay between the start of the ramp gate and the 
start of the plateau gate. 

, t g = time shift in the position of the center of the ramp 

gate with respect to the total two-way delay time of 2h/c. 
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For GEOS-C, the gate width T is matched to the pulsewidth, thus, T_ = 12.5 ns. 

G g 

Since 
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<< i 

Y /h 


and the argument of I (•) in (13) varies as /x» the factor A(t) is almost 
constant over the ranges of integration in (14) and 
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We define the gate weighting function to have unit area, i.e.. 
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T + T + ‘g 


A(T+t ) 
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T - “2 + t g 

The exact shape of the gate weighting function is not known, however, we will 
assume that it is rectangular, i.e., 

1 -T T 
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The time discriminator curve thus reduces to the following form; 
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Equation (16) can be rewritten in a somewhat more convenient form as follows: 
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where S(*) is easily deduced from equation (16). We note from (13) and (16) 
that the A(*) are functions of antenna pointing angle and pattern and altitude 
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while the S(«) are the gate responses to a mean return waveform which is 
free of antenna and pointing angle effects. However, unlike the mean return 
waveform in equation (12), the discriminator output is not a simple product 
of these two factors. 

Equation (16) provides a very rapid means of evaluating the performance 
of split-gate time discriminator systems. It is a significant result in 
itself since it results in a closed form expression. Prior studies of 
discriminator optimization have been hindered by the necessity of accomplishing 
a three-fold numerical integration for E(tg), i.e., a two-dimensional surface 
integration to obtain P r (l) and third integration to determine the gate response 
to P r (x). The above formula has been applied to the Skylab S-193 radar 
altimeter where the pulsewidth is considerably longer (72 ns, 3 dB) and the 
antenna beamwidth smaller (1.5°, 3 dB) . The particular combination of pulse 
width, beamwidth and altitude (435 km) for Skylab results in a mean return 
^hich is neither completely beamwidth nor pulsewidth limited. However, the 
results obtained from (12) and (16) were in excellent agreement with numerical 
integration calculations. We therefore conclude that (12) and (16) have a 
greater range of validity than the approximations would tend to indicate. It 
is interesting to note that from equation (II) it can be shown that the peak 
of the average return power decreases as exp (-4 sin t,f y) £° r very near nadir 
and as l/sin2£ further away from nadir. Such knowledge of the decrease in 
return power as a function of pointing error is very important in designing 
the proper AGC dynamic range of the altimeter receiver. 

4.1.2 Results 

In the following we will present results on time discriminator bias, 
slope and linearity as a function of pointing angle and waveheight. Figure 
4. 3 shows the discriminator curve for the GEOS— C gate configuration (where 

t *■ 

the ramp and plateau gates are separated by 50 ns (T = 62.5 ns)), an rms 
waveheight of 0.1 m, and pointing angles of 0° and 1.5° off-nadir. Apart 
from a shift in the bias point, i.e., the value of t g such that e (t g ) - 0, 
there is no appreciable change in the shape of the curve. We note that the 
linearity range is about -0.3Tg<tg <0.3 Tq Figure 4.4 demonstrates the 
effect of a 1.0 m (rms) surface waveheight on the discriminator curve. 

Comparing these results with those shown in Figure 4.3 indicates a slight 

in bias and a decrease in slope while the linearity range increases 
to -0. 4T t £-0.3T G « The decrease in slope will decrease the equivalent gain 
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of the tracker loop and increase its settling time; however, since the open 
loop bandwidth is unchanged there will be no appreciable increase in tracker 
jitter. Figure 4.5 illustrates the effect of a 2.5 m (rms) waveheight on 
the discriminator curve. Although the linear range of the curve increases, 
there is a marked reduction in slope and the bias is seen to increase both in 
magnitude and sensitivity to pointing angle. Figure 4.6 summarizes the 
resultant altitude bias errors as a function of pointing error and waveheight 
for the GEOS-C Intensive Mode. 

The rather large gate separation in the GEOS-C Intensive Mode tracker 
loop places the plateau gate in that region of the return waveform which is 
sensitive to pointing errors. By reducing this separation, it should be 
possible to reduce the sensitivity of the bias to pointing angle. Figures 
4.7, 4.8 and 4.9 show discriminator curves for a gate separation of 12.5 ns 
and waveheights of 0.1, 1.0 and 2.5 i (rms), respectively. For waveheights 
of 0.1 and 1.0 m, we note a reduction in bias error and sensitivity to pointing 
angle without any significant change in linear range or slope from that obtained 
with the 50 ns gate separation. However, when the waveheight increases to 2.5 m 
(rms), there is a marked increase in altitude bias although the sensitivity to 
pointing angle remains small. In addition, the linear portion of the dis- 
criminator curve is no longer centered about the bias point but has shifted 
to the right of the bias point. The bias errors for a gate separation of 
12.5 ns are summarized in Figure 4.10 and it is noted that the 12.5 ns gate 
separation provides improved performance over the 50 ns separation for low to 
moderate seas but is very poor for high seas. 

Figures 4.11, 4.12, 4.13 and 4.14 show the discriminator curves and bias 
errors for contiguous tracking gates (no separation). We note that the linear 
range of the discriminator curves are" reduced relative to the 50 and 12.5 ns 
gate separation curves. As shown in Figure 4.14 the bias errors for contiguous 
tracking gates are much more sensitive to waveheight than the 12.5 and 50 ns 
configurations . 

4.1.3 Conclusions 

1 in this study we have obtained a concise closed form expression for a 
split-gate time discriminator curve applicable to the GEOS-C Intensive Mode. 
Results for the GEOS-C tracking gate configuration indicate bias errors on the 
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Fig. 4.8 Discriminator Curve for a Gate Separation of 12.5 ns and a g - 1.0 m (H^ ^ - 13.12 ft.) 
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Fig. 4.11 Discriminator Curve for Contiguous Gates and = 0.1 m (H = 1.31 ft.) 



Fig. 4.12 Discriminator Curve for Contiguous Gates and o = 1.0 m (H 1 


13.12 ft.). 
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Fig. 4.14 Altitude Bias Errors due to Waveheight and Pointing Error 
Effects for Contiguous Tracking Gates. 
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order of a few centimeters for the anticipated range of pointing errors and 
waveheights. By reducing the gate separation to one gate width, the bias 
errors may be decreased by about a factor of two for small to moderate wave- 
heights. For large waveheights the bias error resulting from the 12.5 ns gate 
separation increases significantly. In view of this fact*, the optimum tracker 
in the sense of minimum bias error might well be one in which the gate spacing 
is variable and dependent upon the waveheight. For the extreme case of 
contiguous tracking gates, the altitude bias errors are much more sensitive 
to waveheight than the 12.5 and 50 ns gate separations. 

4.2 Pointing Angle Estimation Using the Attitude/ Specular Gate 

In the process of analyzing Sky lab S-193 radar altimeter data, it was 
determined that the pointing angle of the antenna (relative to nadir) could 
be accurately inferred from the shape of the trailing edge portion of the mean 
return. Although the GEOS-C altimeter does not have Sample-and-Hold gates 
located in the trailing edge of the return (as did Skylab), it does have a 
200 ns long integrating gate located in this portion of the return. This 
particular gate has been termed the "Attitude/ Specular" gate; by comparing the 
time averaged output of this gate with the output from the Plateau gate, it 
has been proposed that the pointing angle may be determined. The purpose of 
this section is to investigate how accurately the pointing angle of the antenna 
can be estimated using this technique. 

.Figure 4.1 illustrates how the tracking gates are nominally located 
relative to the idealized return. For both the Intensive and Global Modes, 
the width of the Attitude/Specular gate is 200 ns and the separation between 

it and the ramp-gate is 700 ns. For this study we will make the following 
assumptions; 

(1) tracker jitter, o , , is much~lTess than a pulsewidth and can 

J 

be neglected; 

(2) the tracker gate bias error due to pointing error and sea 
state effects is small relative to the pulsewidth and may 

be neglected, ,i.e. , t = 0; 

8 

(3) the system point target response is Gaussian with a 3 dB 
pulsewidth equal tp 12.5 ns (IM) and 200 ns (GM) ; 

(4) the nominal altitude is 843 km; 

(5) the antenna half-power beamwidth is 2.6°. 
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Since we will be dealing with the Plateau and Attitude/Specular energies, 
the neglect of tracker jitter and bias is certainly justified since these 
effects cause relatively small changes in the location of these gates. 

Let P (t) be the post-video mean return power which has not been 
normalized by the AGC gain. Since we are neglecting jitter, the average 
Plateau gate energy is given by 

Ti + T g / 2 

*p = G AGC f P r (T)dT (18) 

VV 2 


when the GEOS-C altimeter is in acquisition and tracking, the gain of the 
AGC is adjusted so that the average Plateau gate energy is one; thus, 

t i +t g /2 

j P r (x)dT 

VV 2 




The average Attitude/Specular gate energy is similarly given by 


T 2 -T g /2+T 


a/s 


e a/s G AGC 


/'■ 


W 2 


(T)dT 


( 20 ) 


As per GE’s proposed method for determining the pointing angle, the quantity 
of interest is the difference. A, between and e a / s > i*e.. 


( 21 ) 

a/s 

For the Global ttode, T Q - 200 ns and T a/g = 200 ns, so we are integrating over 
comparable time intervals to form e" p and e^ g and A will be less than or equal 
to one. For the Intensive Mode, T G = 12.5 ns and T a / g " 200 ns, and will 
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be roughly 16 (200/12.5) greater than e p ; thus |A|»1. To make A(for both the 
Global and Intensive modes) less than or equal to one, we define as 



( 22 ) 


then using (20) for the Global mode and (22) for the Intensive mode, we see 

that A and A will be on the same numerical scale. It should be noted that 
IM 

if "e . is constrained to a certain numerical range, then there must be some 
a/s __ 

scaling in the hardware in switching from Global to Intensive mode since e^ 
is obtained from different gate widths in the two modes. 

Using the previously derived expressions for the mean return waveform, 
results have been obtained for A and A__. as a function of pointing angle and 
they are shown in Figure 4.15. It is interesting to note that for both modes 
of operation, A is relatively insensitive to pointing angle. The reason that 
the Global mode curve is less sensitive to pointing angle than the Intensive 
mode curve is as follows. For the Intensive mode, the AGC gain is relatively 
independent of pointing angle because the mean return waveform does not vary 
appreciably (over the Plateau gate) with pointing angle. Thus, A.^ determined 
almost completely by the integral over the Attitude/ Specular gate. For the 
Global mode, the AGC gain*decreases with increasing pointing angle while the 
integral over the Attitude/ Specular gate increases. Thus, the product as 
defined by (20) remains essentially constant. 

Since A is the difference of the mean values of two statistical quantities, 
we must know the variance of A in order to state how accurately we can actually 
estimate the pointing angle. There are essentially two error sources involved 
in our estimation of A. The first is due to the noise-like nature of the 
return signal and this is a random error. Unfortunately, it is extremely 
difficult to compute Var(A) because such a computation requires knowledge of 
the autocorrelation function of the non-stationary return waveform process. 

The other important error is due to biases in converting e^ and eT^ to telemetry 
units and then back to- engineering units. In other words, there is an error 
associated with the A/D and D/A conversion process. As per the system speci- 
fication [4], this error is estimated to be ±1% of the recorded value. Thus, 
the bias error for A is given by where 

*The AGC gain and integrated gate values as referred to here are equivalent 
to the quantities defined in (19), (18) and (20) normalized by the factor 
exp[^4sin^C/Y] > (See Sec. 4.1, equations (10) and (11)). 





92 


e. = ±.01 e -(±.01)e , 
b P a/s 


= ±.01-(±.01) (1-A) 


The worst case bias error is seen to be 

E, = ±. 01± .Ol(l-A) 
b 

= ±[ ,02-AJ 

A plot of A±e^ is shown in Figure 4.16. We note that the basic errors involved 

-in converting e" and "e , to telemetry units and back to engineering units 
p a/ s 

limits the angle estimation process to about 0.5 degrees. Thus, we conclude 
that the use of an Aititude/Specular gate to estimate pointing angle is no 
more accurate than about 0.5 degrees. On the other hand, because of the 
relatively large antenna beamwidth, we question whether a 0.5 degree pointing 
error will be an important factor. 
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